Fractal Fidelity as a signature of Quantum Chaos 
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We analyze the fidelity of a quantum simulation and we show that it displays fractal fluctuations 
iff the simulated dynamics is chaotic. This analysis allows us to investigate a given simulated 
dynamics without any prior knowledge. In the case of integrable dynamics, the appearance of 
fidelity fractal fluctuations is a signal of a highly corrupted simulation. We conjecture that fidelity 
fractal fluctuations are a signature of the appearance of quantum chaos. Our analysis can be realized 
already by a few qubit quantum processor. 



Quantum signature of classical chaotic systems have 
been investigated widely in the last decades Within 
the development of quantum information theory, a new 
interest on quantum chaotic systems arose motivated by 
the fact that they are interesting testbeds for quantum 
simulations. Indeed, the quantum simulation of quan- 
tum chaotic systems is possible and some efficient algo- 
rithms, with respect to known classical ones, have been 
found @, Moreover, it has been shown that classi- 

cal simulations of quantum chaotic systems are a difficult 
problem for classical computers due to the high presence 
of entanglement in the fully chaotic regime [5|, |6[ . In a 
more general context the fidelity, the response to a Hamil- 
tonian perturbation of a quantum system, has attracted 
a lot of attention since its introduction by Peres 0]- In- 
deed it has been shown that the fidelity is of fundamental 
importance for the understanding of a system dynamics 
as it has been used to characterize the system intesta- 
bility [8], the feasibility of a quantum simulation H_S 
and of quantum communication protocols [ToL [ill Ml2|], 
the quantum-classical transition [13| , th e sig nature of the 
chaos border [14], the effects of a bath [15] the environ- 
ment induced decoherence [l6[ and the characterization 
of quantum phase transitions [l?], EH, [Hj]. It has also 
been shown that, under given conditions, the fidelity re- 
call classical properties of chaotic systems as its d ecay 
rate is given by the classical Lyapunov exponent [2(il2lf. 

Differently from mentioned previous studies on the 
variance of fidelity fluctuations that characterize the sys- 
tem response (averaged over different initial conditions 
or an ensemble of different perturbations) starting from 
its complete knowledge [26j, we approach the problem 
as for classical complex systems signal analysis. Suppose 
you have a black box you do not have a complete control 
on as, for example, a quantum computer with some given 
but unknown imperfections. Any system outcome, i.e. a 
result of a measurement, might be corrupted: in our ex- 
ample, the result of a quantum computation in presence 
of hardware imperfections. What might be learnt from 
such signal without completely characterize the system 
and/or without affording the cost of repetitions of the 
computation? In this scenario, we focus on the prob- 
lem of finding a clear signature of chaos in a quantum 



systems: it is well known that in classical systems with 
mixed phase space, a fractal region arises at the bor- 
der between integrable islands and the chaotic sea [22^ 
These regions are known to be responsible for fractal con- 
ductance fluctuations as a function of a given param- 
eter (Aharonov-Bohm flux) in open quantum systems 
in the semiclassical regime [23| and it has been shown 
that fractal fluctuations survive also in the deep quan- 
tum regime [241 ] . Moreover, fractal properties of systems 
wave functions have been shown to appear under given 
conditions [25| . Here we show that in a unitary evolution, 
in particular in a quantum computer running a quan- 
tum simulation in the presence of static imperfections, 
fidelity fractal fluctuations naturally arise as a function 
of time. We demonstrate that the fractal dimension of 
such fluctuations strikingly depends on the chaoticity of 
the dynamics: For integrable dynamics the fidelity di- 
mension is integer while for chaotic dynamics it is frac- 
tional. This sensitivity is not restricted to a fully chaotic 
phase space: fidelity fractal fluctuations arise also in the 
chaotic regions of a mixed phase space. We establish this 
connection and we present two possible applications of 
this analisys. First of all it can be used as a testbed 
for a given quantum computer hardware: If running a 
simple (integrable) algorithm fractal fidelity fluctuations 
appears the hardware is not realible as quantum chaos 
has set in [2]. Secondly, we present a method to inves- 
tigate an unknown phase space (Husimi or Wigner func- 
tion) of a general quantum system. Indeed the fractal 
dimension of the fidelity of a given quantum system can 
be used to extract the presence of a chaotic region in the 
system phase space. This is of fundamental importance 
as the simulation of quantum system is one of the most 
general application of quantum computation, however it 
is severely limited by the lack of methods, different from 
the complete wave function tomography, to extract infor- 
mation from the final wave function [27, 28]. In Q, it has 
been shown that the average fidelity can be measured in 
a quantum processor; the sawtooth map algorithm has 
been recently performed on a three qubit NMR quantum 
processor [29j and the fidelity experimentally measured 
on a three qubits quantum processor [30[, thus, the pro- 
posed protocol is at the edge of present day technology. 
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FIG. 1: Fidelity as a function of time for integrable K = — 1 
(red, dashed) and chaotic regime K = y/2 (black, full), with 
e = 10 -4 and n q = 8. Inset: magnification of the behavior 
for long times. 



FIG. 2: Fidelity (Inset) and distribution P(AF) of fidelity 
fluctuations AF = Fi+± — F% for n g = 6, e=10 -3 for chaotic 
(black line) and integrable (red line) dynamics. 



This paper is structured as follows: in Section I we 
present the testbed for our findings, the Sawtooth Map 
and its simulation on a quantum computer in presence 
of imperfections. In Section II we introduce the tools 
to characterize the fidelity fractal fluctuations and we 
analyze the fidelity fractal dimension in different scenar- 
ios. Finally in Section III we introduce the phase space 
tomography via fidelity fractal dimension analysis and 
Section IV is devoted to our conclusions. 



I. SAWTOOTH MAP 

As an example of quantum system with complex dy- 
namics we use the quantum sawtooth map. The cor- 
respondent classical map is characterized by a range of 
very different dynamics depending on system parameters 
varying from completely integrable, semi integrable to 
chaotic dynamics. It is defined by the equations 

{n i+ i = m + k(0i - 7r) 
(1) 
i+ i =0i+Tn i+1 

where (n, 0) are the conjugate action- variables (0 < < 
27r), k is the kick strength and T the time between con- 
secutives kicks. The dynamics is characterized by the 
parameter K = kT: it is chaotic for K > 0, K < —4, 
integrable for K — —3,-2,-1,0 and with mixed phase 
space in the remaining regions [31]. This map belongs to 
the kicked maps family, and it describes a system under- 
going free evolution and subject to a kick every time T. 
The kick strength is proportional to 6 2 . After canonical 
quantization of the action angle variables, the correspon- 
dent quantum map is defined by the Floquet operator 
(time evolution operator for a period T) 

U = e -iTn*^k(§-nf/2^ ( 2 ) 



where n = —id/dO and i/j(9 + 2tt) = ip(6) (we set h = 1). 
The quantization rule is [n, 0] = —i. The map dy- 
namics is still governed by the classical parameter K, 
but a key role is played also by k = K/T which rules 
the quantum-classical transition (k — > oo is the classical 
limit, T = 2tt/N with N = 2 n « the Hubert space size). 
In [2j] an efficient quantum algorithm has been introduced 
to simulate this map and the effects of static imperfec- 
tions, present in any experimental apparatus, have been 
studied [H, 0. The imperfections considered are both 
residual coupling between qubits and fluctuations of the 
single qubits level spacing. As the effects of the static 
imperfections are almost independent from their exact 
expression in terms of operators if there is a stronger lead- 
ing dynamics (the quantum algorithm in our case) here 
we consider only the latter kind of error. The complete 
Hamiltonian that describes the hardware of the quantum 
computer is then 

n q 

H = J2(A + 5 i )al (3) 

i=l 

Here A is the mean qubit level spacing, Si are the fluctu- 
ations of the qubit level spacing taken randomly and uni- 
formly in the interval [— e; e] (constant in time) and erf are 
the Pauli matrices. The Floquet operator in the presence 
of errors is then defined by the quantum gates needed to 
simulate a period of the map ([2]) with the action of the 
Hamiltonian ([3]). The Floquet eigenvalues and eigenvec- 
tors statistics have been studied in [9|, [32[ and it has 
been shown that eventually with growing imperfections 
strength the system ends to be chaotic regardless the 
simulated map dynamics [9[. This crossover, is governed 
by the imperfections strength e and it has been shown 
that in the case of integrable dynamics e c ex riq 5 ^ 2 : this 
threshold for chaos appearance has been estimated from 
the breakdown of perturbation theory and confirmed nu- 
merically [321 ]. 
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FIG. 3: Results of the modified box counting algorithm for 
the fidelity of FigJT] integrable (red triangles) and chaotic 
dynamics (squares). The blue full line is a guide for the eye 
proportional to L 1 ' 36 . Here t* ~ 10 4 . The resulting fractal 
dimensions are D ~ 1.06 for the integrable case and D ~ 1.36 
for the chaotic dynamics. 
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FIG. 4: Fractal dimension of the fidelity fluctuations D as 
a function of Sawtooth Map parameter K, imperfections 
strength from e = 10 -2 to e — 10 -6 (from top to bottom) 
for n q = 8 averaged over Nr < 5 different imperfections real- 
izations. 



II. FRACTAL FIDELITY FLUCTUATIONS 



We study the time fluctuations of the fidelity of a quan- 
tum computation. The fidelity is defined (for pure states) 
as 

F( t ) = \(Mt)\\m)\ 2 , (4) 

that is the overlap as a function of time of the wave func- 
tion computed with the exact time evolution \ijj(t)) and 
the wave function in the presence of static imperfections 
\ip e (t)). The fidelity starts from one, it decays up to a sat- 
uration value and then oscillates around it after a tran- 
sient time t*. 

In Fig. [T] we report two typical fidelities as a func- 
tion of time for a chaotic and an integrable dynamics. 
Notice that the signals are not averaged over different 
static imperfections configuration: they could be the di- 
rect output of an "echo experiment" . We stress the fact 
that the two clear distinct behavior of the fidelity fluctu- 
ations shown in the inset of Fig. [T] reflect the presence, in 
the system classical limit, of a continuum set of typical 
frequencies of the system dynamics in the chaotic case 
differently from the discrete set of harmonics present in 
integrable systems. Even though in Fig. [T] the two curves 
appears pretty different and one could guess the under- 
laying behavior of the system under study, there are case 
where this is not the case. In Fig. [2] we show an exam- 
ple where not only the fidelity time dependence but also 
the distributions of the fluctuations are not distinguish- 
able. In order to extract the information on the under- 
lying dynamics one should use more sophisticated tools 
such as the study of the fractal dimension of the signal 
which detects, as for the classical cases, the presence of a 
"complex" set of typical frequencies of the system in the 
chaotic case and a "regular" spectrum in the integrable 
one. 

The fractal dimension of the signal is measured by 



means of the modified box counting algorithm [33[. In 
the standard box counting algorithm the fractal dimen- 
sion D of the signal is obtained by covering the data with 
a grid of square boxes of size L 2 . The number M(L) of 
boxes needed to cover the curve is recorded as a function 
of the box size L. The (fractal) dimension D of the curve 
is then defined as 

D = — Hrn log L M(L). (5) 

One finds D = 1 for a straight line, while D = 2 for a 
periodic curve. Indeed, for times much larger than the 
period, a periodic curve covers uniformly a rectangular 
region. Any given value of D in between of these integer 
values is a signal of the fractality of the curve. The 
modified algorithm of Ref. [33[ follows the same lines 
but uses rectangular boxes of size L x Aj (A$ is the 
largest excursion of the curve in the region L). Then, the 

number M(L) = l L - is computed. For any curve that 
lays in a plane, a region of box lengths L m ^ n < L < L max 
exists where M oc L D . Outside this region one either 
finds D = 1 or D = 2: The first equality (D = 1) 
holds for L < L m i n and it is due to the coarse grain 
introduced by the discrete map time evolution. The 
second limit (D = 2) is obtained for L > L max and it 
is due to the finite length of the analyzed time series. 
The boundaries L min , L max have to be chosen properly 
for any time series. In our analysis L m i n = T while 
Lmax °° where a increases with n q . In Fig. [3] we 
report an example of this analysis for the long time 
behavior of the two signals of Fig. [TJ We analyzed the 
time series starting from a given time t* to exclude 
the initial transient decay. In Fig. 0] we summarize our 
results reporting the fractal dimension of the fidelity 
for different values of K spanning a wide range of 
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FIG. 5: a,b,c) Fidelity fractal dimension D as a function of 
imperfections strength e averaged over different K for inte- 
gral) le (white) and chaotic regime (black) for n q = 6, 8, 10. 
d) Scaling of the critical imperfections strength e c ~ n~ 2 ' 5 
(white squares, left axis) and limiting fractal dimension D(e — 
0) as a function of n q (black squares, right axis). 



different dynamics (excluding mixed phase space) and 
for different e values. Again, the signature of chaotic 
dynamics is striking. Notice that in the chaotic region 
the fractal dimension of the fidelity fluctuations appears 
to be K independent. Thus, we average over different 
K in both regimes (integrable and chaotic) to study the 
dependence of D as a function of the number of qubits 
and the imperfections strength. In Fig. [5] (a, b, c) we 
show the extended analysis of the fractal dimension of 
the fidelity fluctuations as a function of the number 
of qubits in the quantum hardware and imperfections 
strength for the two distinct dynamics. The difference is 
striking: for chaotic dynamics the dimension is fractal 
for any value of the imperfections strength. On the 
contrary, for integrable dynamics the dimension of the 
fidelity fluctuations is one until a critical value of the 
imperfections strength e c . In Fig. [5] d, we show the 
scaling of e c as a function of the number of qubits in the 
system. Although the scaling is on a very small range of 
parameter, it scales as the critical value for which the 
integrable dynamics became chaotic due to the imper- 
fections effect [32|. It is then clear, that above e c the 
system is chaotic due to the imperfections and regardless 
of the map dynamics: this transition is reflected by the 
appearance of fractal fluctuations. Increasing further 
the imperfections strength eventually the fidelity loses 
all the information regarding the underlying dynamics 
and indeed the fractal dimensions related to the two 
distinct dynamics are equal. The error bars are due to 
the sensitivity of the results to the choice of L min and 
L max . In Fig. [5]d we also show the scaling of the fractal 
dimension D(n) as a function of the number of qubits 
in the limit e — > for the chaotic regime. The signature 




FIG. 6: (Color online) Upper left panel: Husimi function of 
the Sawtooth Map phase space (0, n) for an initial gaussian 
condition inside the chaotic sea (D = 1.34) for e = 2 • 10 -5 , 
n q = 10, t = 2 14 . Other panels: Phase space tomography 
with different resolutions. Fractal dimension of the fidelity 
fluctuations as a function of the initial condition (4 x 4, 8 x 8 
and 16 x 16) for n q = 10 and e = 2 • 10" 5 , Nr = 1. The color 
code goes from blue, yellow to red to represent higher values 
of probability or fractal dimension. 



of quantum chaos is again striking and the difference 
between the two regimes grows with the system size. 

Notice that this analysis can be performed already 
with a few qubits quantum computer (see Fig. [5] a). 
The biggest obstacle to a possible implementation is the 
length of the time series and the correspondent necessary 
long coherent times. However, we checked that already 
with boxes of sizes L ~ 10 in a time series of few hun- 
dred points, the different scalings in the two cases are 
clearly visible. It is also possible, for small e, to per- 
form the same analysis on the initial part of the fidelity 
after subtracting the average decay (data not shown). 
Although this procedure might suffer from errors due to 
fitting procedure of the average behavior, it paves the 
way to an efficient measurement of the fractal dimension 
of the fidelity fluctuations. 



III. PHASE SPACE TOMOGRAPHY 

Finally, we concentrate on the fidelity fluctuations as 
a function of the initial condition in the case of mixed 
phase space. In the upper left panel of Figj6] we report 
the Husimi function of the phase space of the Sawtooth 
Map for K = —2.1, T = 2tt/N where the chaotic sea 
and the integrable islands are clearly visible Q, [34| . The 
initial conditions is a gaussian packet which satisfy the 
minimal indetermination Heisenberg condition centered 
inside the chaotic sea. In the lower panels of Figj6] we 
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show the result of the phase space tomography with the 
proposed method with a discretization of 4 x 4, 8 x 8 and 
16 x 16 initial conditions. For every initial condition, a 
gaussian packet centered in (no,#o), we computed the 
fractal dimension of the fidelity fluctuations D(no,#o) 
and plotted the resulting contour plot from D = 1 
(black) to D = D max (red). We found again that the 
dimension of the fidelity depends on the underlying 
dynamics: a fractional dimension is found for the chaotic 
sea while it is integer when the dynamics lies inside 
the integrable islands. The result is a partition of the 
phase space in chaotic and integrable regions which 
can be refined increasing the number of different initial 
conditions (notice that the resolution of the Husimi is 
of the order of N 2 = 2 20 points). Thus, this analysis 
can be used as a tool to scan an unknown phase space 
of a given Hamiltonian dynamics to discern between 
integrable and chaotic phase space regions. 



IV. CONCLUSIONS 

In conclusion, the fidelity time fluctuations are a signa- 
ture of quantum chaotic dynamics. They reflect, as for 



classically chaotic signals, the presence of an "almost- 
continuum" of frequencies in the system dynamics com- 
pared to the "discrete" set of typical frequencies of in- 
tegrable systems. A similar behavior has been already 
observed in the context of quantum communication in 
spin chains [lOj. We propose this kind of analysis to in- 
vestigate an unknown phase space of a complex quantum 
system and to test the effects of static imperfections of 
quantum hardware. We stress that the sensitivity of the 
fidelity fluctuations of a simulated system dynamics is 
valid also in the limit of imperfections strength going to 
zero, i.e. they will be present in any experimental quan- 
tum hardware in case of chaotic or mixed phase space dy- 
namics. Similar results are likely to be found in presence 
of time dependent imperfections (errors) in the quantum 
system along the lines of [35[ . 
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